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We study exchange coupling in Si double quantum dots, which have been proposed as suitable 
candidates for spin qubits due to their long spin coherence times. We discuss in detail two alternative 
schemes which have been proposed for implementing spin qubits in quantum dots. One scheme 
uses spin states in a single dot, and the interdot exchange coupling controls interactions between 
unbiased dots. The other scheme employs the singlet and triplet states of a biased double dot as 
the two-level system making up the qubit and exchange controls the energy splitting of the levels. 
Exchange in these two configurations depends differently on system parameters. Our work relies on 
the Heitler-London approximation and the Hund-Mulliken molecular orbital method. The results 
we obtain enable us to investigate the sensitivity of the system to background charge fluctuations 
and determine the conditions under which optimal spots, at which the influence of the charge noise 
is minimized, may exist in Si double quantum dot structures. 



I. INTRODUCTION 

The use of quantum-mechanical two-level systems to 
construct quantum bits, or qubits, for the storage and 
manipulation of information has attracted considerable 
attention in the past decade. Computational algorithms 
relying on the principles of quantum mechanics such as 
unitary evolution, entanglement, and coherent superpo- 
sition are expected to offer novel solutions to problems 
such as factoring, searching, and simulating quantum- 
mechanical systems. 1 Among the many physical systems 
that have been proposed as candidates for qubits, solid 
state systems are particularly promising because of then- 
potential for scalability, and spin-based qubits in solid 
state systems are regarded as ideal candidates due to 
their observed long coherence timesi^ 

A common essential feature of the most prominent 
spin-based quantum computing (QC) architectures is the 
employment of exchange interaction between localized 
electrons. For example, in the original Loss-DiVincenzo 
(LDV) proposal 4 - exchange interaction is pulsed on and 
off by electrical control and produces two-qubit gates 
such as SWAP and controlled-NOT gates. The Kane 
proposal^ of using nuclear spins of 31 P in Si as qubits 
employs donor electron exchange to mediate nuclear spin 
interactions and produce two-qubit gates. In the more 
recent singlet-triplet (ST) qubit proposal,— controlled 
exchange splitting between singlet and triplet states leads 
to single-qubit gates. Exchange interaction is electro- 
static in nature so that it can be strong and can produce 
fast gates (on a time scale as short as tens of picosec- 
onds), and it can be controlled electrically by chang- 
ing voltages applied to the gates defining the quantum 
dots. In other words, exchange gates are easy to control 
and allow straightforward interface with existing micro- 



electronic devices. A potentially important shortcom- 
ing of an exchange-based QC architecture is also related 
to its electrostatic character: turning on the exchange 
interaction could potentially make the system vulnera- 
ble to electrical fluctuations in the environment (back- 
ground charge noise j 9 - - — gate noise ^ electron-phonon 
interaction ; 13 i 14 etc.), which can lead to qubit decoher- 
ence. Therefore, it is imperative that the exchange inter- 
action be characterized carefully, so that its magnitude is 
known in a physical system, and possible ranges of con- 
trol parameters for which the influence of the external 
charge noise is minimized (so-called optimal points) can 
be identified,! -!^ 

One of the most promising semiconducting materials 
to host a spin quantum-information processor is silicon, 16 
which has outstanding spin coherence properties in the 
bulk. 17 Natural silicon contains only a small number of 
nuclear spins (5% of 29 Si, which has spin 1/2), and can 
be further isotopically purified. The strength of the hy- 
perfine interaction between a conduction electron (or a 
donor-bound electron) with each individual nucleus is 
also weak; 1 ^ As such the hyperfine-interaction induced 
electron-spin decoherence, dominant in III-V materials 
such as GaAsj 1 ^— 7 is suppressed in Si. Furthermore, 
Si has a small spin-orbit interaction 2 ^ and no piezoelec- 
tric interaction, which together means slower electron 
spin relaxation due to phonons. 29 Silicon does have a 
drawback in its somewhat complicated conduction band 
structure with six equivalent minima. This orbital de- 
generacy could lead to difficulties such as the atomic- 
scale position dependence of exchange interaction be- 
tween donor electrons^ - At present there are several pro- 
posed Si quantum-computer architectures, including ar- 
chitectures based on donor electron or nuclear spins in 
a Si:P system^ ' 31 ' 32 single electron spins in gate-defined 
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quantum dots in Si/SiGe^ or Si/Si02(Ref. 33) quantum 
dotSj2^~— and SiGe nanowires^i 

The development of Si single-electron devices suitable 
for quantum-information processing is a relatively recent 
phenomenon because of the difficulty in identifying and 
growing appropriately lattice-matched barrier materials 
(as compared to Gai-^Al^As for GaAs), which help min- 
imize carrier scattering near the interface of a hetero- 
junction, and difficulties in fabricating and characterizing 
single-donor devices at the atomic length scale. Never- 
theless, exciting experimental progress has been made in 
both donor-based devices^— and gate-defined quantum 
dots^-^.-il 

Motivated by the experimental attempts to realize spin 
qubits in a Si mctal-oxide-semiconductor field-effect tran- 
sistor (MOSFET) j^ir— here we consider the prospects 
for creating spin qubits in a lateral double quantum dot 
(DQD) structure fabricated in a silicon MOSFET by 
lithographic patterning. We calculate the exchange cou- 
pling and tunnel coupling of electrons in DQD structures 
in a Si/Si02 system under the assumption of a reason- 
ably large valley splitting. We first analyze the appli- 
cability of the Heitler-London (HL) and Hund-Mulliken 
(HM) models for calculating exchange splitting in a Si 
double dot, then perform our calculations for experimen- 
tally achievable sizes of dots. We study two types of 
DQD configurations. One is a symmetric DQD, where 
the single-electron levels of the two dots are on resonance 
and are appropriate for the exchange gates of the orig- 
inal LDV proposal. The other type is a biased DQD, 
where the two dots are voltage-biased to the vicinity of 
the resonance between the two-dot singlet and one of the 
doubly occupied singlet states. This is a configuration 
appropriate for the ST qubit architecture. 

Compared to the calculation of exchange coupling in 
GaAs^ ' 12 i 15 ' 48 " — there are some distinctive features for 
such a calculation in Si. First of all, in a Si MOSFET 
structure the Coulomb interaction is enhanced by the 
proximity of a lower- k dielectric, while the kinetic en- 
ergy is reduced by the larger effective mass. Therefore 
correlation effects should be much stronger in a Si dou- 
ble dot as compared to a GaAs structure, shrinking the 
range of validity for the simpler models of exchange cal- 
culation (as we will discuss in Sec. Ill CJ) . Indeed, iden- 
tifying these ranges is one of the primary aims of the 
present calculations. Another important distinction is 
the presence of conduction-band valleys in Si. Bulk Si 
has six degenerate conduction-band minima. Confine- 
ment and uniaxial strain in the z direction split these 
six valleys into a two-dimensional(2D) manifold of lower 
energy and a four-dimensional manifold of higher energy 
that is separated from the lower one by several tens of 
meV(Ref. 31 and 51) and is typically neglected in dis- 
cussions of the lowest-energy states. Furthermore, the 
interface potential typically splits the two lower-energy 
valleys by the energy on the order of a fraction of an 
meVj£2r— The presence of the valley degree of freedom 
complicates the energy spectrum and spin structure of 



the two-electron states in silicon quantum dots^ 5 - — and 
can potentially affect the effectiveness of Si as a host for 
spin-based quantum-information processing. 

In this paper we consider situations in which the val- 
ley splitting A (between the ground and the first valley 
eigenstates) is reasonably large so that consistent loading 
of the nondegenerate ground valley-eigenstate in a quan- 
tum dot or a double dot can be achieved. We also confine 
ourselves to the situation where the valley composition 
of the ground valley eigenstate is uniform throughout the 
DQD. 57 Under these conditions we can perform our cal- 
culation of exchange within a single valley because the 
two-electron Hamiltonian Eq. ([T]) does not have any in- 
tervalley matrix elements: all the single-particle terms 
are slowly varying in space, and all the Coulomb matrix 
elements (except one exchange integral) vanish exponen- 
tially ~ exp(— L 2 /afy where L is the ^-direction confine- 
ment and is on the order of 3 nm while ao is the Si lattice 
constant of 5.43A, so that the exponential factor is ex- 
tremely small. The lone exchange integral that does not 
vanish exponentially is very small in any case (< 0.2/ieV 
for reasonable QD sizes) so that we take it as zero as 
well. Since the double-dot Hamiltonian does not lead to 
any coupling to the higher-energy valley states even if an 
excited valley is nearby, it will not produce any modi- 
fication in the single-valley exchange we calculate, and 
the basic physics we consider here is the same as in a 
single- valley system such as GaAs. 

In this study we also investigate the sensitivity of 
exchange-coupled spin qubits to environmental electri- 
cal fluctuations. As we mentioned above, finite-exchange 
coupling makes spin qubits susceptible to electrical noises 
from sources such as gate electrodes and background 
charge fluctuations. Background charge fluctuations 
arise due to trapped electrons in defect sites at semicon- 
ductor interfaces. Movement of nearby trapped charges 
affects quantum dot systems in two ways: rise or fall 
of the barrier between the dots and fluctuations in the 
two confining potentials of the dots. Due to the pos- 
sible importance of charge noise (Si02, the most com- 
monly employed barrier material for silicon hetero junc- 
tions, is a well-known source of charge noise), the condi- 
tions for the existence of the so-called optimal points (in 
analogy to the optimal point studied for superconduct- 
ing charge qubits, where charge-noise-induced dephas- 
ing is minimized 5 ^), where the exchange splitting for the 
two-electron states is the least sensitive to environmental 
charge fluctuations, need to be investigated. In this work 
we therefore also calculate the sensitivity of the qubit to 
the charge noise and identify the optimal points in the 
qubit parameter space in which the qubit is less sensitive 
to fluctuations in gate voltagesi 1 ^ 11 ' 15 ^ 9 

The goal of this paper is a qualitative (and at best 
semi-quantitative) description of the exchange-gate be- 
havior in Si DQD qubit architecture as a function of 
the relevant physical parameters: dot separation, con- 
finement potential, dot size, interdot bias, and external 
magnetic field. Since the details of the actual Si DQD 
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confinement in real samples are unknown and likely to be 
extremely complicated, we use a simple physically moti- 
vated model potential to study the problem. Given the 
simplicity of our model, our calculation of exchange also 
employs well-known approximation schemes (i.e. Heitler- 
London and Hund-Mulliken) which should be qualita- 
tively valid in various limits. At this stage of the devel- 
opment of Si spin qubits, it is not particularly meaning- 
ful to carry out a fully numerical Schrodinger-Poisson- 
configuration-interaction modeling of the exchange-gate 
energetics, because the details of the structures are sim- 
ply not known and will probably remain unknown for 
some time at the desired level of accuracy to justify a 
completely computational first principles approach. Our 
simpler approach with relative computational ease en- 
ables us to provide a thorough exploration of parame- 
ter dependence of the Si DQD exchange-gate energetics, 
which should be helpful in this early stage of experimen- 
tal development. In this sense, our work should be con- 
sidered as the Si equivalent to the corresponding GaAs 
qubit theories carried out in Refs. andlla 

The outline of this paper is as follows. Section II in- 
troduces the formalism that will be used throughout the 
paper, which includes the model potentials for double 
quantum dots, the procedures followed in the determina- 
tion of exchange coupling, and the discussion of their lim- 
its of validity. In Sec. Ill the two architectures employed 
to make spin qubits in quantum dots arc discussed. Sec- 
tion IV contains our results on exchange coupling and 
sensitivity to noise. Section V contains a summary and 
conclusions. 



II. THEORETICAL FORMALISM 

A. Hamiltonian and double-dot-model potentials 

Consider a DQD at the Si/Si02 interface with growth 
direction along the £-axis. The effective-mass Hamilto- 
nian for two electrons in a DQD is then 



V(ri) is the confinement potential chosen to approxi- 
mate the real potential for two dots separated by a po- 
tential barrier, E is the electric field applied in the x- 
direction (the interdot axis direction of the DQD), and 
the last term is the electron Zeeman energy (here we have 
assumed that both electrons are experiencing the same 
magnetic field), where for Si we can take g C R = 2. 

In the above hi we have omitted the spin-orbit cou- 
pling terms 6 -^ since these are expected to be very small 
in silico n 28 ! 61 compared to GaAs. The spin-orbit terms 
in general induce an anisotropic (Dzyaloshinskii-Moriya- 
type) exchange interaction between the two spins ; 60 ' 62- — 
which could affect the performance of two-qubit exchange 
gates. 66 While an experiment aiming at detecting this 
kind of anisotropic exchange interaction was proposed, 66 
these effects have not yet been clearly shown to be rele- 
vant to experiments on GaAs (where they are expected 
to be much stronger than in Si). Thus, at this stage of 
research on Si DQDs we choose to neglect them in this 
work. 

The potential V(ri) describes the DQD system of two 
dots located in the xy-plane which are tightly confined in 
the z direction. The exact form of the confinement along 
the z direction is irrelevant here as long as the char- 
acteristic localization length in the ^-direction is much 
smaller than in the x and y directions. This is fulfilled 
in Si/Si02 DQD where the z-confinement length is on 
the order of 3 nm, while the confinement length in the 
x and y-directions is at least several times larger. The 
strain at the Si/Si02 interface and the confinement in 
the £-direction lifts the sixfold valley degeneracy in bulk 
silicon and leaves the ground valleys doubly degenerate. 

This ground-valley degeneracy is further split by an 
amount A due to the interface potentialr&^^i^i 6 -^ 
which allows us to focus only on the lowest-energy valley 
eigenstate. At this stage we take the 2D approximation 
in which the finite confinement length in the ^-direction 
is neglected^ The reliability of this approximation will be 
ascertained in Sec. Ill CI We will use three different mod- 
els of the DQD potential. The first is a quartic model^ 
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V(r) = ^muj, 
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(x 2 - d 2 ) 2 + y 2 



(4) 



where £ = 1,2 labels the two electrons, hi is the single- 
particle Hamiltonian, k — (est + Esio 2 )/2 is the effec- 
tive dielectric constant including the image charge in the 
adjacent Si02, and r%2 is the distance between the two 
electrons. The single-electron Hamiltonian is 



hi 
f, 



% + V{Ti) + eExi + g eS fi B BS iz , (2) 

-L[ Pi -f A (r,)] 2 , (3) 
2m c 



where m = 0.191m e is the transverse (isotropic in the xy 
plane) effective mass for electrons in the z and — z val- 
leys (the valleys contributing to the lowest-energy valley 
eigenstate) in Si, A = B(—y, x, 0)/2 is the vector poten- 
tial for magnetic field B applied along the ^-direction, 



while the second is a biquadratic model^ - 

V(r) = ^mcu 2 (Min[{(x-d) 2 ,(x + d) 2 }] + y 2 ) . (5) 

Note that both of these potentials provide infinite con- 
finement (i.e. we have V(r) — > oo as r goes to infinity). 
A third, more realistic potential, will be introduced in 
Sec. Ill C[ where we will use the comparison of calcula- 
tions with all three model potentials to gauge the limits 
of validity of our theoretical approach. 

For the quartic potential from Eq. ((4]) the two dots 
are separated by a distance 2d and the central barrier 
height is V& = |towqgP . For the biquadratic potential 
from Eq. the distance between the dots is the same, 
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while the central barrier height is Vb = ^rauS^d 2 , four 
times larger than for the quartic potential. Note that 
for these potentials the barrier height and the interdot 
distance are not independent variables. On the other 
hand, in experiments there are typically three indepen- 
dent voltages controlling the DQD potential (two gates 
controlling the depth of the left and right dot potentials, 
and an additional gate separately controlling the barrier) . 
However, most often only two variables (the left-/right- 
dot gate voltage) are changed simultaneously, which we 
model with a nonzero in-plane electric field E in Eq. ([1} . 

Each of the model potentials is designed in such a way 
that around the centers of the two dots at x = ±d (for 
E = 0) the confinement is approximately parabolic so 
that the ground state single-particle wave-functions are 
harmonic oscillator states in the isolated dot limit. More 
precisely, near the centers of the two dots we have 

V(x « ±d) « \m<4 \{x T d? + y 2 } = V R/L (r) , (6) 

and the approximate single-dot lowest-energy eigenstates 
( < 5 L/fl( r ) are solutions of the equation 



[Ti + V L/R (r)](p L/R (r) = hw (p L/R {r) 



(7) 



with <Pl/r{v) — (r\L/R) being the Fock-Darwin states 
centered on the right/left dot^ 



fL/R 



exp[=F ^ ] exp[ 



b((x±d) 2 +y 2 ) 
2a 2 1 

where the first exponential is the magnetic phase for a 
displaced orbital, a = ^h/mojQ is the Fock-Darwin ra- 
dius, and b= y/l + u l /lu 2 is a measure of the degree of 
magnetic confinement, where lol = eB/2mc is the elec- 
tron Larmor frequency. The overlap of the right and left 
single-particle orbitals is 



I = (L\R) = exp 







[-( 




a 2 \ 





(9) 



For example, at zero magnetic field the overlap is I — 
cxp(-d 2 /a 2 ). 

A ST qubit is operated in the parameter regime where 
the lowest-energy (1,1) singlet (one electron in each dot) 
is close in energy to the lowest energy (2,0) singlet state 
(both electrons in the left dot, as shown in Fig. [IJ. For 
small quantum dots the on-site Coulomb repulsion is 
quite strong, so that the DQD needs to be strongly biased 
and asymmetric. We model such a biased DQD by adding 
an in-plane electric field E = Ex along the inter-dot axis 
to the Hamiltonian. While adding E displaces the two 
minima of the potential corresponding to the centers of 
the two dots for both model potentials, other responses 
are significantly different for these models. Specifically, 
when an electric field is added to the quartic potential, 
the potential minima of V(x) are altered by E and thus 
the curvatures of V(x) near the potential minima are also 
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FIG. 1: (color online) The potential of an asymmetric double 
quantum dot in the ^-direction (hujo = 6 meV). The solid line 
corresponds to the biquadratic potential while the dot-dashed 
line corresponds to the quartic potential. The ratio of the 
detuning to the single-particle energy near the anti-crossing 
(explained in the text) is e/ftwo ~ 4. 



changed, which in turn changes the single-dot confine- 
ment energies of the L/R states (for example, in Fig. Q] 
we show a case in which there is no minimum on the 
right side, only a saddle point). In addition, as shown in 
Fig. [TJ the height of the barrier between the two dots is 
strongly affected by E in the quartic model. In short the 
quartic model of the DQD becomes rather cumbersome 
when the two dots are biased. These observations lead 
us to adopt the biquadratic model for the calculations 
with finite bias field E. It is reasonable to expect that 
the behavior of this model is closer to the experimen- 
tal situation, in which the primary effect of changing the 
gate voltages is to introduce an offset between the ground 
state energies of the two dots. Changes in the barrier 
height (i.e., the tunnel coupling between the dots) and 
in the confinement energy in each dot are secondary ef- 
fects. 

Using the biquadratic potential from Eq. ([5]) with an 
added electric field E, the two minima of the potential 
(corresponding to the centers of the two dots) are shifted 
by the same distance /a, with / given by 



/ 



eEa 
hcua 



(10) 



The single-dot orbitals are now given by Eq. (JS) with d 
replaced by d — fa for the right dot, and — d replaced by 
—d — fa for the left dot. 

In experiments on asymmetric DQDs the commonly 
used parameter is the so-called detuning energy e. It 
parameterizes the changes of the DQD potential along 
a specific line in the experimentally obtained charge- 
stability diagram^ of the DQD. A change in e corresponds 
to a change in the difference between the single-particle 
ground states in the two dots. Using the biquadratic 
potential model, we get e = 2ffkuo(d/a) — 2eEd. 



5 



B. Theoretical approaches for calculating the 
exchange splitting 

In a uniform magnetic field, if two electrons interact 
only through the Coulomb interaction, the spin eigen- 
states are the singlet (|5) = A=\ tl — it)-, w hh to- 
tal spin 5 = 0) and triplet (|T ,+,_) = 4j| U + it), 
| ff), and | II)), which have total spin 5 = 1) states. 
The corresponding spatial parts of the wave functions are 
symmetric and antisymmetric under exchange of the two 
electrons, respectively. For quantum computing, where 
orbital degrees of freedom are frozen out, the energy split- 
ting J = Et — Es between the unpolarized triplet state 
| To) and the ground singlet state is the key physical quan- 
tity and we will refer to it as the exchange splitting or, 
in general, exchange interaction. 

To calculate the exchange splitting J, the two-electron 
Hamiltonian should be diagonalized in the singlet and 
triplet subspaces, respectively (with symmetric and an- 
tisymmetric orbital-state bases) , which amounts to a full 
configuration interaction (CI) calculation. Since there is 
an infinite number of two-electron states, all practical CI 
calculations are truncated. Depending on the degree of 
truncation, there is a hierarchy of models for the calcu- 
lation of exchange interaction in a two-electron double 
quantum dot. 

Interestingly, the reliability of a CI calculation is not 
completely determined by the size of the finite basis left 
after the truncation^ 

One particular test of validity of the given approach 
to calculating J is the Lieb-Mattis theorem^ which dic- 
tates that at zero magnetic field the singlet must always 
be the ground state so that J > 0. The condition for ap- 
plicability of this theorem in 2D and three-dimensional 
(3D) cases is for the confinement potential to be sep- 
arately symmetric, i.e., symmetric under exchange of 
each of the x, y, and z coordinates of the two particles. 
This is the case here, since we use model confinement 
potentials which are separable. At finite B either sign of 
J is P ossibleAi2^a 



left and right orbitals given in Eq. ([9]). Heitler-London 
wave functions are essentially educated guesses of the 
two-electron ground states. They are reasonable wave 
functions for quantum dot states at large interdot sepa- 
rations. 

In the HL basis the exchange splitting is the energy dif- 
ference between the unpolarized To triplet and the singlet 
states 



Jhl = (T \H\T ) - (S\H\S) , 
which can be written as 



2l 2 



l-/ 4 



w v 



Dq - j^E 



(12) 



(13) 



where W v is the single- particle contribution due to the 
confinement potential 



W v = ( L(1)J2(2) 



1 



L(l)R(2)--L(2)R(l)j 



(14) 



withw(l,2) = F(l) + V r (2)-Vi ? (2)-V L (l). The analytical 
expressions for W v for quartic and biquadratic potentials 
are given in Appendix. Physically, W v represents the 
gain in kinetic energy for the singlet state arising from 
the fact that electrons in this state are spread out more 
evenly across the two dots. Dq is the direct Coulomb 
interaction contribution 

A, = (L(1)R(2)\C\L(1)R(2)) , (15) 

and Eq is the exchange integral 

E = {L(1)R(2)\C\L(2)R(1)). (16) 

Both Coulomb terms are positive. Together they favor 
the triplet state since in an antisymmetric state the two 
electrons avoid each other to minimize their interaction. 
Both Coulomb integrals Dq and Eq can be done analyti- 
cally in 2D (Ref. 9) using the states given in Eq. ©. In 
units of the effective Rydberg Ry* = me 4 /2H 2 k 2 (approx- 
imately 44.76 meV in a Si/SiC>2 system) we obtain 



Da 



bd 2 



-bd 2 /a 2 



(17) 



1. Hettler- London approximation 

The simplest approach for calculating the exchange 
splitting J in a DQD is the HL method^ ' 12 i 70 ' 71 where 
the ground two-electron double-dot singlet and triplet 
states are built from single-dot single-electron states in 
the simplest possible combination. The spatial parts of 
the normalized singlet and triplet states are then given 
by 



|5/T) = 



\L(l)R(2))±\L(2)R(l)) 



(11) 



where \L/R(i)) are the single-dot ground state orbitals 
for the ith electron, given by Eq. (|5J| for the symmet- 
ric DQD case, and 1= (L\R) is the overlap between the 



Eq = l 2 



(6- 



6- 



(18) 



where Iq is the zeroth-order modified Bessel function, a = 
a/og and d= d/a^ are the Fock-Darwin radius and the 
half interdot distance expressed in units of the effective 
Bohr radius = h 2 K,/me 2 (approximately 2.11 nm in 
a Si/SiC>2 system). Since d/d = a/d, we can drop the 
tildes when the ratio appears in Eqs. (JTTJ) and (TT51) . The 
contributions of Dq and Eq to J can be of comparable 
magnitude. For example, at large distances (d^>a), we 
have Eq ~ l 2 D d/a at B = 0, and Eq « ID at high B fields 
(when b»l). The I /I 2 factor in front of E a in Eq. {T3J) 
makes the negative contribution of the exchange integral 
Eq significant, and possibly larger in absolute magnitude 
than the direct term Dq (this clearly occurs at high B, 
when J can become negative). 
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2. Hund-Mulliken Approximation 

The HL approximation neglects both the contribu- 
tion of the higher-energy orbitals from the two dots and 
the possibility of double electron occupation on either 
of the dots. The latter restriction is removed in the 
Hund-Mulliken (HM) method^^^ISJl in which the 
doubly-occupied states (built out of the ground-state 
single-dot orbitals) are included in the basis set. This 
is of particular importance when we consider the asym- 
metric case of a biased DQD, which is relevant for the 
experiments on singlet-triplet qubits. In that system the 
mixing between the delocalized singlet 5(1, 1) and one 
of the doubly-occupied states [e.g. 5(2,0)] by the tun- 
neling of an electron between the two dots is the main 
mechanism that brings about the singlet-triplet splitting 
J [(nL,n R ) denotes the number of electrons in each dot]. 

In the HM method we take into account the doubly- 
occupied (2, 0) and (0, 2) states. We still employ only 
the lowest-energy single-dot orbitals, so that the doubly- 
occupied states must be spin singlets (their spatial parts 
being symmetric). These two doubly-occupied states, to- 
gether with the (1, 1) spin-singlet state, form the basis of 
three singlet states in the HM method, which have an 
anti-symmetric spin part but a symmetric spatial wave 
function. Conversely, the spin-triplet state, with a sym- 
metric spin part, has only one possible spatial configura- 
tion since the electrons cannot occupy any higher-energy 
orbital states. It is convenient to use the basis of or- 
thonormalized single-dot states, which we define follow- 
ing Ref. [3 as <5> L/R = {iPl/r - g(p R /L)/y/l - 2lg + g 2 , 
where g = (1 — y/l — l 2 )/l. Then the spatial wave func- 
tions of the three singlets and the (1,1) triplet are given 
by 



*i /fi (ri,r 2 ) = $L/K(ri)$L/jt(r a ) 
(l.i)/ 

S/T ' 



(19) 



1 



(ri,r 2 ) = -^=^ L (r 1 )$H(r 2 )±$fl(ri)$i(r 2 )j(20) 

where the superscript d denotes the doubly occupied 
states. The Hamiltonian written in the 5(2,0), 5(0,2), 
5(1, 1), T(l, 1) basis takes the form 



H = e R + e L 



X 
s[2t 



X y/2t \ 

U + e V2t 

y/2t V s 

V T J 



(21) 



where 



£R 



x = (*t /R \c\K/L) u = (n, R \c\n, R ) (22) 

V s = (9 s \d\9g) V T = {*t\C\Vt) 
t = ($ L/R \h\$ R/L ) + j=(* s \C\* d L/R ) =t' + w 



In the above, s l and e R are the single-particle energies in 
the left and right dots, e is the detuning parameter, X is 
similar to a interdot Coulomb exchange integral, U is the 
on-site (Hubbard) Coulomb repulsion, Vs and Vr are the 
Coulomb energies for the (1,1) singlet and triplet states, 
with ^s/t built out of orthogonalized single-electron or- 
bitals, and t is the interdot tunneling matrix element 
modified by the Coulomb matrix element w (f is the 
bare tunneling). All these can be expressed in terms 
of the analogous matrix elements between the original, 
nonorthogonal, left and right orbitals fL/R (these bare- 
orbital matrix elements are termed Do, Eq, Uq etc). Ex- 
plicit formulas are given in Appendix. 

Diagonalizing the Hamiltonian given in Eq. (|2ip can 
then give the exchange interaction which is the energy 
difference between the two lowest eigenvalues [one cor- 
responds to a singlet, the other the (1,1) triplet]. For 
the symmetric case in zero B field, the lowest eigenvalue 
£s_ is dominated by the energy of the 5(1, 1) spin singlet 
state while modified by the two doubly occupied states. 
On the other hand, for the asymmetric (biased) DQD, 
es_ is dominated by a mixture of the 5(2, 0) and 5(1, 1) 
singlet states [we assumed here that the bias is such that 
5(2, 0) state is the lower-energy doubly occupied state, 
i.e., the potential of the left dot is lowered to attract an 
additional electron]. 5(0,2) state can then be dropped 
from our consideration. 

We define the critical bias detuning e c as the bias at 
which the 5(1,1) state is on resonance with the dou- 
bly occupied 5(2,0) state when tunneling t is neglected. 
Then from Eq. (pi]), the energy of 5(1, 1) is 



£ s(i,i) — £ R + £ l + Vs , 
and the energy of the 5(2, 0) state is 
£s(2,o) — 2el + U . 



(23) 



(24) 



Therefore the critical bias detuning e c corresponds to the 
energy difference e c = U — Vs, at which £s(i,i) = £ s(2,o)- 
When e <C e c , the ground state is composed mostly of 
5(1, 1); when e ^> e c the 5(2, 0) state becomes the ground 
state. Around the anticrossing region where e ~ e c the 
(1, 1) singlet and (2, 0) singlet states are strongly mixed. 

We note here that the bias regime in which we are in- 
terested in our HM calculations is where the S(l,l) sin- 
glet is the dominant component of the ground singlet 
state. This is the regime where single-qubit gates are 
performed for an ST qubit. The high-bias regime, where 
S(2,0) is the dominant component of the ground state, 
is important for initialization of an ST qubit, but is not 
considered in this study. 



3. Beyond the Hund-Mulliken model 

Beyond the Hund-Mulliken model one can build CI cal- 
culations with various degrees of truncation. For exam- 
ple, by including the single-electron p orbitals, one can 
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introduce anisotropy into the electronic states^ which 
allows the electrons to spread their wave functions better 
in order to minimize the total energy (therefore a better 
chemical bond), as compared to the HL and HM mod- 
els, which are based solely on s orbitals. The inclusion 
of higher orbital states also allows a better account of 
the correlations between the two electrons, leading again 
to lower total energy. On the other hand, it is also im- 
portant to note again that a finite-basis CI calculation 
is not necessarily always better than the simple HL and 
HM calculations.- - Furthermore, the simple and analyt- 
ical HL and HM models have the advantage of providing 
transparent physical pictures. Therefore in the present 
calculation we use only the HL and HM models, with 
particular emphasis on identifying the range of validity of 
these models under various circumstances. We note here 
that a CI calculation was recently done to explore the 
highly biased regime where the doubly occupied 5(2, 0) 
is the ground stated This is a complementary study to 
what we explore in this work, in analogy to the case of a 
helium atom versus a (biased) hydrogen molecule. 



C. Reliability of exchange calculations 

The HL and HM methods used in this paper gener- 
ally underestimate the electron correlations in the dou- 
ble dot, and results obtained using these approaches are 
more appropriately thought of as order-of-magnitude es- 
timates. However, these methods do capture the basic 
physics of exchange coupling in DQD systems and calcu- 
lations performed under these approximations can serve 
as guidelines for initial development of silicon-based DQD 
structures. In this section we discuss the expected regime 
of applicability of these methods. 

We first clarify the relationship between the valley 
splitting A and the other energy scales in the problem as 
we work within a single-valley approximation. To have 
reliable control over the initialization process, when the 
electrons tunnel into the quantum dot from an outside 
reservoir, it is necessary for A to exceed the thermal 
broadening of the reservoir Fermi level: A ksT. If 
Coulomb interaction can couple different valleys, it could 
pose additional requirements on A. Indeed, Coulomb 
interaction can connect charge densities from different 
valleys, whether on the same dot or on different dots. 
However, Coulomb matrix elements involving overlaps of 
states from different valleys are strongly suppressed^ due 
to the large wave vectors k z = ±0.85 (27r/ao) that sepa- 
rate the valleys, where ao = 5.43 A is the lattice con- 
stant of silicon. Therefore, as long as the two-electron 
system is prepared in the ground valley eigenstate, nei- 
ther changes in quantum dot potential (which is slowly 
varying in space) nor Coulomb interaction can couple 
the electrons to the higher-energy valley states. In other 
words, the only requirement on A in our calculation is 
A> fcflT for initialization. 

One of the main sources of problems in HL and HM 



calculations is the fact that the single-dot states used 
to build the two-electron bases can be unreliable, espe- 
cially for small interdot distances and/or for small barrier 
heights. This is related to the well-known failure of the 
HL method when the overlap I between the L and R 
states becomes too large. It was shown 9 that when the 
parameter c= e 2 / (^afiwo) = \f^j2d (proportional 

to the ratio of the single-dot Coulomb energy to the con- 
finement energy) becomes larger than 2.8 in 2D, the HL 
method gives the unphysical result of J < at zero B 
for small d. This occurs also in 3D case with isotropic 
mass^i only at higher values of c > 5.8. Thus, above a 
certain size for the quantum dots (as measured in terms of 
a and determined by the ratio between Coulomb and con- 
finement energy) , the HL and HM approaches will surely 
fail qualitatively below a threshold in d. Compared to a 
GaAs double dot with the same d and a, in a Si/Si02 
system the Coulomb energy is larger due to the reduced 
screening caused by the adjacent Si02 layer while the ki- 
netic energy is smaller due to the larger effective mass. 
Both these changes lead to smaller <Zg . Consequently, the 
critical value of c is reached in smaller dots in Si/Si02 
(larger hujQ to compensate for smaller n, from a different 
perspective) compared to those in GaAs, reducing the 
range of dot sizes that can be accurately modeled by the 
HL and HM methods. 



Clearly, the case of closely-spaced dots can only be 
addressed by more sophisticated quantum-chemical ap- 
proaches where a larger basis of single-particle orbitals is 
usedi 12 i 50 ' 73 Here we choose to perform HL and HM cal- 
culations for weakly coupled dots. The criterion used to 
ascertain the qualitative applicability of these simple ap- 
proaches is for the exchange splitting J to increase mono- 
tonically with decreasing d at zero magnetic field. We 
define the critical distance 2d c as the interdot separation 
at which the sign of the derivative dJ/dd changes, and J 
starts to decrease toward zero and negative values at even 
smaller d. In Fig. HJa) one can see that the HL and HM 
methods fail at values of d c which are in the range of pos- 
sible interdot separations in experiments. The critical d c 
as a function of confinement energy is shown in Fig.[^b), 
from which we can see that for interdot distance of 2d~ 50 
nm we can only consider dots with hujQ > 5 meV, which 
corresponds to Fock-Darwin radii a<9 nm for Si/Si02- 



While for d > d c the exchange splitting changes ac- 
cording to the expected trend, we still need to establish 
the reliability of HL and HM calculations quantitatively. 
For this purpose we have adopted an approach proposed 
in Ref. [74, in which the reliability of the HL approach 
was significantly improved by employing a more realistic 
DQD potential and choosing an appropriate variational 
form of the single-electron orbitals for this potential. The 
modified potential is illustrated in the inset of FigJ5Ja). 
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Analytically it is given by 

|mw 2 (Min {(x - d) 2 , (x + d) 2 } + y 2 - r 2 ) , 



MV-HL 



V 



d) 2 ,(x + d) 2 } 



y 2 < r 2 



if Min {(a; 
0, if Min {(x - d) 2 , (x + d) 2 } + y 2 > r 



(25) 



where r is the radius at which the confinement poten- 
tials on both dots are truncated, different from the Fock- 
Darwin radius of the dot. In our calculations we choose 
d/r in the range of 0.48-1.6. The depth of the dot po- 
tential that we choose in the calculations below allows 
at least three single-particle energy levels to exist in 
each dot. Fig. HJa) shows that the potential function in 
Eq. (|2l)|) yields exact harmonic confinement around the 
quantum dot centers and becomes zero at larger distances 
from the centers. The single-particle wave functions have 
to be modified in order to reflect the fact that the poten- 
tial barrier around the double dot is not infinite. We 
consider single-particle wave functions that are centered 
at the minimum of each dot to take the following varia- 
tional form: 



<p(p) 



Ai exp 



A 2 exp 



2a? 



ap 
2ai 



, P < P 



, P > P 



where p denotes the distance from the center of the single 
dot. This piece- wise wave function involves new param- 
eters A\, A2, /i, a and /3. However, the normalization 
condition and the boundary conditions (continuity of the 
wave function and its derivative) at (x =p d) 2 + y 2 = p 2 
would leave us two independent parameters, which we 
choose to be f3 and p, the first of them related to the ra- 
dius of the electron wave function in the dot and the sec- 
ond being the radius at which the change of decay behav- 
ior of the wave function occurs. The parameters j3 and 
/i are obtained variationally by minimizing the single- 
particle energy in a single dot, the potential of which is 
given by 



V(P) 



\mu 2 {p 2 - p 2 ), p<r 
0, p > r 



This matched variational (MV) wave function adopted 
here gives a larger overlap compared to the Gaussian or- 
bitals used previously since the wave function has an ex- 
ponential tail. The variational procedure also allows the 
MV orbitals to be better matched to the potentials of 
the two dots. Furthermore, the total potential at large 
interdot distances is now simply a sum of the two po- 
tentials corresponding to the two dots, in analogy to the 
potentials of two atoms that make up a molecule, which 
suggests that the single-dot MV wave functions are a bet- 
ter starting point for HL and HM calculations. 

The same approach can be used in the case of nonzero 
B field. For a single dot we use the gauge in which 
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FIG. 2: (a) Exchange energy as a function of interdot dis- 
tance at confinement energy 8 meV (c w 4.18, a w 7.1 nm). 
The solid line corresponds to the Heitler-London calculation 
with the MV wave functions as in Ref. [fU (r « 24 nm and 
f3 « 0.089). The dotted line is HM calculation in 2D GaAs 
system with the same size of dot. The long-dashed line is a 
3D calculation using the HL method with confinement radius 
3 nm in the z-direction. The short-dashed line corresponds to 
the 2D HM method with the quartic potential model while the 
dot-dashed one is the 2D HM approach with the biquadratic 
potential model. The inset shows the modified double quan- 
tum dot potential V from Eq. 1251 (b) Critical interdot dis- 
tance (2d c ) as a function of confinement energy using the 
quartic potential model. The solid line corresponds to the HM 
method while the dashed line corresponds to the HL method. 



A = B(—y,x,0)/2, which produces an effective parabolic 
potential ^muj 2 L p 2 and a term proportional to the L z , 
the projection of the angular momentum on the z axis. 
Since we use a trial wave function <p(p) given above which 
depends only on p, only the B-induced harmonic confine- 
ment needs to be taken into account. The calculation is 
analogous to the previous one. 

The improvement of the results in the MV approach 
over the HL calculation with Gaussian orbitals obtained 
in Ref. [74] is consistent with the results from Ref. [Hoi , 
where HL results for a biquadratic DQD potential were 
compared with a multiorbital calculation. The MV ap- 
proach was also shown to be in qualitative agreement 
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with the analytically interpolated expressions for J in 
the double-donor potential ) 74 ' 75 where one can use the 
exact results for J obtained using the Gorkov-Pitaevskii 
methodi 7 - 6 - 

These results indicate that the MV calculation should 
be a good approximation to the exact result for exchange 
splitting, at least for larger d. In Fig. EJa) we show 
that using the MV approach we can calculate J for more 
closely spaced dots than in the case of HL or HM approx- 
imation with Gaussian orbitals. The dip in J(d) calcu- 
lated by the MV method at d ~ 26 nm is most probably 
an artifact due to using the HL method with the confine- 
ment potential having a kink between the dots — a simi- 
lar feature was seen in HL calculation with a biquadratic 
potential in GaAs^ Furthermore, the values of J ob- 
tained using the different methods and potentials exhibit 
the same behavior at d > d c . The spread in values of 
J obtained with different potentials and wave functions 
is about an order of magnitude, but since the J itself 
changes by orders of magnitude with quite small rela- 
tive change in d, this feature does not affect the qualita- 
tive information which we can gain from our calculations. 
The main point here is that at d > d c various methods 
of calculating J agree, strengthening our assertion than 
at these interdot distances the results of HL and HM 
approaches are semi-quantitatively reliable. In Sec. IIVI 
we will present the results of calculations for symmet- 
ric and asymmetric (biased) DQDs using these methods, 
with the quartic potential and the modified MV poten- 
tial for a symmetric DQD, and the biquadratic potential 
for the asymmetric DQD. We do not use the biquadratic 
potential for the symmetric DQD because the unrealistic 
kink in the central barrier overly suppresses the interdot 
overlap and biases the result in favor of the triplet state, 
resulting in the failure of the HM calculation at much 
larger value of d c . 

In Fig. [U[a) we also show the result of the HL calcu- 
lation with 3D Gaussian orbitals, i.e., the wave func- 
tions from Eq. ([8) multiplied by the Gaussian envelope 
exp(— z 2 /2a 2 L ) in the z direction. Here a± is the Fock- 
Darwin radius in the z direction, which in our calcula- 
tion is a± = 3 nm. As expected from the lowering of 
the Coulomb energy in the 3D calculation due to the 
increased spread of the wave functions, the critical inter- 
dot distance d c is smaller in the 3D case, thus increas- 
ing somewhat the range of validity of the HL and HM 
approaches* 7 ^ However, since the confinement in the z- 
direction is generally stronger than in the xy-plane and 
the Fock-Darwin radius associated with it smaller, mak- 
ing the spatial integrals anisotropic, we cannot obtain a 
closed-form expression for the Coulomb contribution to 
the exchange interaction. 

For comparison we have also performed the calcula- 
tions for GaAs using the same interdot distance and the 
same Fock-Darwin radii a as the corresponding Si dots 
(an example is shown in Fig. [5]). The resulting exchange 
splittings are typically an order of magnitude larger than 
the ones obtained in Si. However, since decreasing d by 



about 10% allows one to gain an order of magnitude in 
J, there is no practical obstacle in obtaining the same 
values of J in Si and GaAs. 

In summary, in this section we have shown that the HL 
and HM methods are valid in smaller parameter regimes 
in a Si/Si02 system as compared to GaAs DQDs due 
to the enhanced Coulomb interaction and larger effective 
mass in Si/Si02- These are facts one needs to be aware 
of when calculating electronic states in Si. One conse- 
quence for experiment is that the QDs in Si have to be 
made smaller and closer in order to have large exchange 
splittings. From comparison of calculations with differ- 
ent methods and different model potentials we draw the 
conclusion that the HL and HM calculations are semi- 
quantitatively accurate at moderately large interdot dis- 
tances. 



III. SPIN QUBIT ARCHITECTURES IN 
QUANTUM DOTS 

Our study of electron exchange in a DQD is mainly 
motivated by two approaches to spin quantum compu- 
tation using quantum dots. In the original spin QC 
proposal, which we refer to as Loss-DiVincenzo (LDV) 
architecture; 4 - an array of quantum dots is used as a set 
of spin qubits. Each dot is singly occupied and the spin 
of the electron in each dot constitutes the qubit. Single- 
qubit operations can be implemented by rotating the sin- 
gle spins using a magnetic field. Two-qubit operations, 
which are needed in order to implement logical gates, can 
be realized by means of a pulsed (gate-controlled) ex- 
change interaction between neighboring dots by lowering 
the inter-dot potential barrier and allowing the electron 
wave functions to overlap. Remote two-qubit operations 
can be realized by swapping the two qubits next to each 
other, again assisted by pulsing on and off exchange in- 
teraction between neighboring quantum dots. 

The LDV scheme involves a system of symmetric dots 
since, in general, no electrical bias is applied between 
the adjacent dots at any time. Indeed small bias about 
the symmetric point will only reduce the exchange split- 
ting J ~ The model potential is symmetric and the best 
choice for the two-dot scenario is the quartic potential 
introduced in Sec. |TTJ 

The simplicity and elegance of the LDV architecture 
have motivated extensive experimental studies of sin- 
gle spin properties such as coherent control, coherence, 
and measurement^' 8 - ' 37 ' 46 ' 77- — The experimental prob- 
lems encountered in these explorations, such as diffi- 
culties in reliable and fast spin initialization and single 
spin ESR control, have led to theoretical proposals where 
pure electrical control of encoded spin states can be used 
for universal quantum computation. 82,83 In particular, 
a series of experiments ^ 79 ' 84 ' 85 have shown the viability 
of two-spin encoding using spin singlet and unpolarized 
triplet states. 

For a singlet-triplet qubit, the exchange splitting J 
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plays the same role as the Zeeman splitting does for a 
single-spin qubit. Turning it on produces a phase differ- 
ence between the singlet and triplet states so that dif- 
ferent superposition states can be realized. For the ex- 
perimentally realized ST qubit the double dot is in 
a strongly biased regime so that (1,1) and (2,0) singlet 
states are close to resonance in the absence of tunneling. 
The exchange splitting is turned on by starting deep in 
the (1,1) regime [where (1,1) singlet and triplet states are 
nearly degenerate ground states of the DQD, indicating 
that the interdot tunnel coupling is small] and increasing 
the mixing between the (1,1) and (2,0) singlet states due 
to tunneling In other words, the exchange splitting 
here is dominated not by the Coulomb interaction but 
by tunneling-induced mixing between the levels detuned 
by e. 

One of the potential drawbacks of working with a 
highly biased double dot is that charge noise affects the 
bias linearly, which in turn leads to fluctuations in the ex- 
change splitting J.— This amounts to a dephasing mech- 
anism for the singlet-triplet qubit if J is finite. There- 
fore in the current study we also investigate whether it is 
possible to reduce the sensitivity of this system to charge 
fluctuations. 



IV. RESULTS 

Based on the above discussions of the reliability of the 
HL and HM approaches within the specific regime we 
are considering, we calculate the exchange splitting be- 
tween the triplet and singlet states in a coupled DQD 
near a Si/Si02 interface. Our calculations are done for 
experimentally achievable sizes of dots. We investigate 
the dependence of exchange on the confinement energy, 
interdot distance, and the magnetic field in the symmet- 
ric (Loss-DiVincenzo) case, and the sensitivity of the ex- 
change coupling of the two LDV qubit to charge noise. 
For the asymmetric case (singlet-triplet qubit), we de- 
termine the dependence of the exchange coupling on the 
interdot detuning within the Hund-Mulliken model, and 
we discuss the presence of the optimal points with respect 
to different kinds of charge noise. 



A. Loss-DiVincenzo architecture 

Our discussion of the Loss-DiVincenzo architecture re- 
lies on the quartic model potential for reasons discussed 
in Sec. Ill CI where we showed that this potential gives 
physical results for J in a wider range of interdot dis- 
tances compared to the biquadratic model. 

In Fig. [3] we plot J as a function of ^ at zero mag- 
netic field. The results here are straightforward to un- 
derstand. As the interdot separation increases, the over- 
lap I — exp(— d 2 /a 2 ) decreases as a Gaussian. Corre- 
spondingly the exchange coupling J decreases according 
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FIG. 3: (color online). Exchange coupling as a function of 

,2 

=y using the quartic model potential. The solid and dashed 
lines correspond to HM and HL with hwo = 8 meV (a ~ 7.1 
nm) respectively while the dot-dashed line corresponds to HM 
with huJo = W meV (a«6.3 nm). 



to Eq. (fTB")) , and the decrease is dominated by the I 2 pref- 
actor for J. Fig. [3] shows that with a confinement energy 
of 8-10 meV, exchange can reach 1/^eV in the range of 
interdot distance of about 30 nm. In experiment, the 
confinement energy will typically be smaller than 8 meV, 
which implies that the actual exchange coupling at the 
Si/Si02 interface could be orders of magnitude larger. 
We note that J ~ 1 /xeV is already well within energy 
resolution of current generation experiments,-. It cor- 
responds to a gating time of about a nanosecond. The 
coherence time of a single spin in GaAs quantum dot has 
been measured in a spin echo experiment to be about 300 
ns at low magnetic field (B<0.1 T), 3 and longer coher- 
ence times are expected at higher B fields i 20 ' 21 i 23 ~ — Con- 
sequently, even a J of 1 /LteV magnitude should be enough 
for quantum-computation purposes within the LDV and 
ST architectures. 

Figure 0] shows J as a function of confinement en- 
ergy hujQ. Increasing the confinement energy changes 
the central barrier height (14 = imw^d 2 ) between the 
two dots, but more importantly it squeezes the Fock- 
Darwin radius since a = ^/h/rnujo, which has a similar 
effect on the exchange coupling as varying the interdot 
distance. When the confinement energy is larger than 
10 meV (with a < 6.5 nm), the exchange coupling be- 
comes extremely small at interdot distances above 30 
nm since I 2 now is smaller than 10~ 10 . On the other 
hand, when the confinement energy is smaller than 6 
meV (with a > 8 nm), we obtain an unphysical triplet 
ground state for interdot distances smaller than 32 nm. 
The reason for this unphysical behavior, which has also 
been observed for interatomic calculations before, lies in 
the failure of HL and HM models to properly account for 
electron correlations^ 

In Fig. [S] we plot J as a function of a uniform applied 
magnetic field B in the z direction. The DQD system 
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FIG. 4: (color online). Exchange coupling as a function of 
confinement energy using the quartic model potential for two 
dots separated by 40 nm. The solid line corresponds to the 
HM approach while the dot-dashed line corresponds to the 
HL method. 



shows a clear transition from an antiferromagnetic (sin- 
glet) ground state at low field to a ferromagnetic (triplet) 
ground state at high field. This trend, which appears here 
in (single- valley) Si DQDs, is consistent with previous re- 
sults obtained for GaAs DQDs.^^^ In essence, as the 
external field increases, the interdot overlap decreases so 
that the overall magnitude of J goes down. However, 
the quantum interference brought about by the magnetic 
phases in Eq. © causes the magnitude of the exchange 
Coulomb term to increase relative to the direct Coulomb 
term (as shown in the inset of Fig. [5] for the HL and MV- 
HL calculations, where Eq/1 2 — Dq increases monotoni- 
cally with the applied magnetic field. The singlet-triplet 
crossing occurs at the field where the quantity plotted 
in the inset reaches the value of -1), so that at larger 
fields the triplet configuration is favored energetically, 
as shown in Fig. [S] More involved numerical calcula- 
tions performed for GaAs DQD have shown that further 
changes in the two-electron ground state are possible at 
higher fields. 48 However, whether the triplet will remain 
the ground state at even higher magnetic fields cannot 
be determined by our simple model calculation — indeed 
the inset of Fig. [5] seems to indicate that there will be 
no further singlet-triplet crossings at higher fields. The 
exchange coupling will go to zero asymptotically with in- 
creasing B, since the interdot overlap becomes so small at 
large B that the two dots are effectively fully separated. 

Figure [B] shows J as a function of the interdot distance 
d at finite B. The top panel shows the results obtained 
using the quartic model potential and HM method, while 
the bottom panel shows the results obtained using the 
MV approach. As shown in Fig. for a given d the ex- 
change coupling J changes its sign with increasing B at 
specific interdot distances. Changing d while keeping B 
fixed again changes the magnetic phase factor, leading 
to possibly nonmonotonic dependence of J on d. Inter- 
estingly, the behavior of J as a function of d is different 
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FIG. 5: (color online). Exchange coupling as a function of 
magnetic field at confinement energy 8 meV and interdot dis- 
tance 30.24 nm. The solid line corresponds to the MV-HL 
calculation, and the dashed (dot-dashed) lines are HM (HL) 
calculations with the quartic potential. The inset presents 
the value of (Do — Eo/l 2 )W v [see Eq. (|13[) ] as a function 
of magnetic field, which indicates the relative importance of 
Coulomb interaction compared to single-particle kinetic and 
potential energy terms. 



in these two models. When the MV approach is used, 
the magnitude of J is monotonically decreasing with in- 
creasing dot separation, with the exception of a small 
range of B around 6.5T (shown in the inset). Within 
this range the curve of J versus d is at first positive, 
then drops below zero and goes through a minimum. For 
these few values of B there is thus a point on the curve 
where dJ/dd = 0. On the other hand, when the quartic 
model potential is used, J is negative at smaller d and 
relatively low B, and it has a maximum over a wide range 
of B. Since the negative J at smaller d in the quartic po- 
tential calculation can actually be due to the failure of 
the HM method at small interdot distances, and in view 
of the wider range of applicability of the MV method 
(see Fig. [2]), we expect that the trend shown by the J(d) 
curve produced by the MV method as the magnetic field 
increases is more likely to match realistic experimental 
situations, while caution must be exercised in interpret- 
ing results obtained using the quartic model potential. 

When two electron spins are exchange coupled, the 
spins are exposed to the electrostatic environment. Fluc- 
tuations in the electrostatic potential will lead to fluc- 
tuations in the exchange splitting J, which in turn 
leads to gate errors ^ In general we can write A J — 
^2,ApJ/dVi) AVi, where the partial differentiations re- 
flect the sensitivity of the double dot system to the par- 
ticular fluctuation, and the Vi are different sources that 
contribute to the barrier or the interdot bias.— For ex- 
ample, dJ/dVb, with Vb being the barrier in the model 
potential, should tell us how much voltage fluctuations 
in gate electrodes that define the DQD could affect the 
exchange, and it is important to study dJ/dVb and ways 
to minimize it. 
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FIG. 6: (color online). Exchange coupling as a function of 
interdot distance at confinement energy 8 meV with different 
magnetic field. Fig. (a) shows the result of quartic model. 
Fig. (b) shows the result of MV model(r ~ 24 nm) The inset 
shows the result of exchange coupling as a function of interdot 
distance at B = 6.5 T. 



While the interdot distance in the model quartic po- 
tential is directly related to the central barrier height 
Vb = mui^d 2 /8, such a simple relation generally cannot 
be established for more realistic potentials. It certainly 
does not hold for the truncated potential used in the MV 
calculation, the results of which should be more reliable. 
We would like to point out that the derivative dJ/dd 
quantifies the response of J to fluctuations of the exter- 
nal potential which predominantly cause the change in 
the interdot distance, i.e., they are even with respect to 
yz plane. The inset of Fig. [3] indicates that there is an 
optimal point (at which dJ/dd = Q) with respect to such 
fluctuations predicted by the MV calculation at B = 6.5 
T. The corresponding plot of dJ/dd is shown in Fig. [7] 
Outside the small range of magnetic field around 6.5 T, 
the magnitude \dJ/dd\ is simply monotonically decreas- 
ing with increasing d. We note that dJ/dd in the absence 
of a magnetic field shows a qualitatively similar behavior 
as those with B ^ 6.5 T and there is no optimal point in 
that case. 
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FIG. 7: (color online), d J/dd (the derivative of J with respect 
to the interdot distance) as a function of interdot distance at 
a confinement energy of 8 meV and at B = 6.5 T, using the 
MV approach within the HL method. 



B. Singlet-triplet qubit 

The description of the ST qubit necessarily requires 
taking into account the doubly occupied state in one dot, 
since tunneling between (1,1) and (2,0) singlet states is 
the physical process that lowers the energy of the singlet 
relative to the triplet state, as we discussed in Sec. Mil 
Therefore this problem can only be tackled using the HM 
approximation. In modeling a biased DQD we make use 
of the biquadratic potential with an electric field added 
along the x direction. This choice is due to the fact that 
the curvature of the quartic potential is altered consid- 
erably by the addition of the electric field, and the min- 
ima at the centers of the two dots have rather different 
shapes, as can be seen in Fig. [TJ On the other hand, the 
problems with the biquadratic potential in the symmetric 
case (i.e. the failure to obtain J>0 occuring at interdot 
distances larger than when using the quartic potential) 
are not as severe in the biased case: the physics which 
controls J in the regime of parameters relevant for the 
operation of the ST qubit is not qualitively sensitive to 
the details of the potential shape. 

The exchange coupling in a biased DQD is the energy 
difference between the two lowest eigenvalues of the HM 
Hamiltonian in Eq. (|2"Tj) . In Fig. [5] we show the three 
lowest eigenvalues of the HM Hamiltonian as a function 
of the detuning energy near the (1, l)-singlet to (2,0)- 
singlet anticrossing. The evolution of the DQD spectrum 
near this charge transition is parametrized by the detun- 
ing e = 2(d/a)fhuo = 2eEd [see Sec. II. A and Eq. (10)], 
which is a linear function of both the interdot distance 
and the electric field. For relevant experimental manipu- 
lations it is natural to assume that the change in electric 
field is the main driving force behind changes in detun- 
ing. However, one could also keep f^E fixed, and vary 
d in order to vary the detuning. As e increases in Fig. [51 
the ground state of the DQD changes gradually from the 
(1,1) spin-singlet state to the (2,0) spin-singlet state. 
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FIG. 8: (color online). Energy spectrum as a function of 
bias near the anticrossing. The confinement energy is 6 meV 
and E=0.6 V//xm. We have used the asymmetric biquadratic 
potential. 



The anticrossing point at which e = e c can be viewed, 
approximately, as the point where the composition of 
the ground state changes from predominantly (1,1) to 
mostly (2,0). We focus on the spectrum near this an- 
ticrossing, specifically for e < e c , which is the range of 
detuning used in qubit manipulations [the regime where 
5(2, 0) is the ground state is used only for initialization 
of the ST qubit]. For these values of e we calculate the 
exchange coupling, tunnel coupling, and we investigate 
the sensitivity to charge noise of the asymmetric DQD. 
In the following we use a confinement energy of 6 meV, 
where c = 4.84, the Fock-Darwin radius a = 8.16 nm and 
d/a > 2. The exchange coupling in the biased DQD can 
be then tuned in a wide range (0 ~ lQO/ieV) near the 
anticrossing point. 

Fig. [9] shows the exchange splitting J and tunnel cou- 
pling t as functions of e near the anticrossing point. The 
exchange splitting depends on both tunnel coupling and 
detuning, and it becomes very large near the anticross- 
ing (resonant) point between the 5(1, 1) singlet and the 
5(2, 0) states. In the upper (lower) panel we show the 
results when detuning is changed via variation of applied 
electric field (interdot distance) . Comparison of the two 
panels shows that the electric field and the interdot dis- 
tance have quite different effects on J and t. As the 
electric field increases (upper panel), the tunnel coupling 
increases slightly (less than 5% in the given range of de- 
tuning) because the barrier height relative to one of the 
dots decreases slightly. In addition, the overlap between 
neighboring orbitals does not depend on the electric field, 
so that the Coulomb terms do not change at all, and t is 
only affected by slightly different barrier height. There- 
fore, the change in J in this case (varying by a factor 
of 7) is dominated by the change in the detuning, since 
J ^2t 2 /(e c — e) in the regime where S(l,l) is the domi- 
nant component of the ground state, away from the anti- 
crossing point. Conversely, while increasing the interdot 
distance leads to a larger e as well (lower panel), it also 
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FIG. 9: (color online). Exchange coupling (solid line) and 
tunnel coupling (dashed line) versus the bias near the (1,1)- 
(2,0)-singlet resonant point. The upper panel (a) presents the 
results when change of bias is by varying the electric field, 
with the interdot distance fixed at 40 nm and confinement 
energy fttJo = 6 meV. The lower panel (b) presents the results 
when change of bias is by varying the interdot distance d with 
the electric field fixed at E=0.6 V / fim and confinement energy 
huo — 6 meV. 



lowers the overlap between the left and right orbitals and 
exponentially reduces the tunnel coupling. The exchange 
coupling in this case is determined by the competition 
between the energy level detuning and the tunnel cou- 
pling. The curve of J as a function of bias is no longer 
monotonous but decreases first and then rises as shown in 
the lower panel of Fig. |H1 The presence of this extremum 
in J(d) dependence leads to the existence of an optimal 
point which we discuss below. 

Figure QJJ shows J as a function of B at different de- 
tunings. The overall decrease in J is again due to the 
magnetic squeezing of the orbital states and suppression 
of interdot overlap, as in the symmetric dot case. There is 
also a transition from antiferromagnetic (singlet) to fer- 
romagnetic (triplet) spin coupling for the ground state, 
yet this feature is far less dramatic than in the symmet- 
ric case. The expanded view in the inset shows that the 
exchange splitting switches sign from positive to nega- 
tive when the magnetic field exceeds 9 T although the 
negative J is very small in magnitude. A possible reason 
for this suppressed singlet-triplet crossing is that in a bi- 
ased DQD the exchange is dominated by the variations 
in detuning and tunnel coupling. While Coulomb inter- 
action modifies these factors, it is not as important as in 
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FIG. 10: (color online) Exchange coupling as a function of 
magnetic field for two values of the bias detuning e with con- 
finement energy of 6 meV and interdot distance of 40nm. The 
dashed line corresponds to e = 22 meV while the solid line 
corresponds to e = 23 meV. The inset gives an enlarged ver- 
sion of J(B) at large magnetic field. 

a symmetric DQD. Therefore, the effect of competition 
between direct and exchange Coulomb integrals is also 
less prominent. 

In an asymmetric DQD the effect of background charge 
noise on J is governed by the dJ/dd and OJ/dE deriva- 
tives. The first of them parameterizes the influence of 
fluctuations in the DQD potential that are even with re- 
spect to yz plane, while the second controls the influence 
of the odd changes. In Fig. [IT] dJ/dd is plotted at rela- 
tively low and relatively high magnetic fields as a func- 
tion of bias e (the interdot distance 2d is varied). Panel 
(a) shows d J/ dd versus the bias detuning at small mag- 
netic fields. There is a sweet spot on this graph at which 
dJ /dd = 0, which is also seen in the absence of a mag- 
netic fields This is a point where exchange splitting has 
a local minimum and the reason is the competition be- 
tween a decreasing tunnel coupling and a detuning that 
is approaching (l,l)-(2,0) resonant point which we dis- 
cussed before (see also the lower panel of Fig. [SJ . On the 
other hand, in panel (b) we show the results at higher B 
fields. At Bk8 T a second optimal point appears. The 
exchange splitting has its local maximum at this point, 
which is due to the contribution of magnetic phase. If 
the magnetic field is larger than 10 T, dJ/dd is always 
positive and has its minimum value near the anticrossing. 

In Fig. [T2]we show dJ/dE as a function of detuning 
at an interdot distance of 40 nm. The sensitivity to 
charge noise, encapsulated by dJ/dE, increases mono- 
tonically as the bias increases, and there is no optimal 
point. The physical reason for this monotonic depen- 
dence is straightforward: as bias increases, the charge- 
distribution difference between the singlet and triplet 
states acquires a stronger dipolar characteristic [in the 
(1,1) configuration this difference is electric quadrupo- 
lar], so that the DQD system becomes more susceptible 
to environmental charge fluctuations. 59 
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FIG. 11: (color online) d.J/dd as a function of the bias e by 
varying interdot distance, with a fixed electric field of E = 0.6 
V/fim and a single-dot confinement energy of huuo = GmeV. 
Panel (a) shows one sweet point dJ /dd = at small magnetic 
field(_B < 6 T). While Panel (b) shows that there are two 
sweet spots around B — 8 T but no sweet point for B > 10 
T. 

We note here that while clear separation of physical 
effect that changing E and d have on J is a feature spe- 
cific to the model biquadratic potential, in which the ap- 
plied electric field does not change the interdot distance, 
our result highlights the fact that the dependence of J 
on the potential fluctuations can exhibit local extrema if 
the relevant potential fluctuation influences the effective 
interdot distance. 



V. SUMMARY AND CONCLUSIONS 

We have calculated exchange coupling of two elec- 
trons in double quantum dot structures in a Si/Si02 
heterostructure within the single-valley approximation. 
This approximation is valid as long as the valley split- 
ting A 3> ksT, which is necessary for reliable initializa- 
tion of electron-spin states for the purpose of quantum 
computation. Having analyzed the reliability of Hcitlcr- 
London and Hund-Mullikcn models for calculating ex- 
change splitting in a Si double dot, we investigated sym- 
metric and asymmetric double quantum dots of experi- 
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the possibility, as was done for GaAs DQD structures^ 
that the individual dots in the Si DQD structure may 
have more than one electron each to determine whether 
the limit of one-electron per dot is essential in Si qubit 
architectures or any small odd number of electrons per 
dot could work under experimental conditions. Finally, 
the question of what happens if the valley splitting is not 
large (i.e. compared to the temperature) is an important 
issue as discussed recently in Ref. [13, and must also be 
dealt with in a fully quantitative manner. 



FIG. 12: (color online). d.J/dE versus the bias for different 
magnetic fields. The bias is varied by altering the electric field 
E. The solid line corresponds to B = 1 T while the dashed 
one corresponds to B — 3 T. 
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mentally achievable sizes, with our focus on two-electron 
exchange splitting and its sensitivity to charge noise. 
Throughout this work we have also used a matched vari- 
ational approach^ to extend the range of applicability of 
our results. 

For the Loss-DiVincenzo architecture we find that, for 
the same interdot separation as in GaAs, exchange is gen- 
erally smaller because of the larger effective mass and the 
resulting smaller wave-function overlap. In the presence 
of a magnetic field, we find that for different confinement 
potentials the trend exhibited by J as a function of d as 
the magnetic field increases, and thus the existence and 
location of sweet spots, is different. Modeling symmetric 
double quantum dots in a magnetic field with HL and 
HM methods requires a thoughtful choice of model po- 
tential so as to avoid unphysical results. We find that 
optimal point (at which the sensitivity to charge noise is 
reduced) exists for J as a function of d for a small range 
of magnetic fields. 

For singlet-triplet qubits, we find that the behavior of 
the exchange as a function of electrical bias and of mag- 
netic field is qualitatively similar to GaAs. We identify 
an optimal point at B = 0. We find that the sensitivity to 
the noise depends on whether the potential fluctuations 
are changing the interdot distance d or not, i.e. whether 
they are odd or even with respect to the yz-planes. 

Future work ought to include a Poisson-Schrodinger 
modeling of the actual Si DQD structures, yet this must 
await the successful fabrication of Si DQD structures so 
that more details of the gate-induced electrostatic con- 
finement associated with the lithographic processing be- 
come known quantitatively. It will also be useful to go be- 
yond our simple Heitler-London and Hund-Mulliken the- 
ories of the exchange coupling using quantum-chemical 
configuration interaction corrections (Ref. [73| is an ex- 
ample along that line) although we do not expect any 
qualitative modifications of our results with more sophis- 
ticated numerical computations. One might also consider 
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Appendix A: Matrix elements in the Hund-Mulliken 
approach 



In the following, we list the matrix elements in Eq. (|21[) 
and Eq. (|22[) in units in which the energy is measured in 
hujQ and length is measured in Fock-Darwin radius a — 
\J h/mujQ. The single-dot harmonic well-ground eigen- 
states ifiL m for dots centered at ±d are given by Eq. |[8j . 
The dimensionless parameters are then: b = \/l + uj^/uJq 
(related to relative strength of magnetic confinement ver- 
sus harmonic confinement), where u>l is the Larmor fre- 
quency ojl = eB/2mc; c = yjir /2e 2 / (nafoujQ) (the ra- 
tio of Coulomb to confinement energy); and d = d/a 
(the ratio of half interdot distance to Fock-Darwin ra- 
dius). In addition, the overlap is given by I = (L\R) = 
exp [dp (1/6 — 2b)] ; and the parameter associated with 
the applied electric field is f — eEa/hujQ. 

In order to convert the expressions below into the units 
of Ry* and a|j one has to multiply all the matrix elements 
by huo = 2/a 2 (note also that c = \Jit/2 a). 

The matrix elements in Eq. (|2"Tj) are as follows: 
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we have the on-site Coulomb repulsion term: 



d 



± 



( e -6(<io+/) 2 (_ 1 + e 4Mo/)) 



2(1 — Z 2 ) 

T(do + /) evi\Vb(d + /)! ± (do - /) erf\Vb(d - f) 



Ua = (L(1)L(2)\C\L(1)L(2)) 
= (R(1)R(2)\C\R(1)R(2)) 
= cVb , 



(A7) 



to which we can relate the interdot Coulomb exchange 
interaction 



1 



, 2 f 2 d e- b ^+^ 2 



■\l 2 [-2b+f 



b + di- J —- 



2 

2d e- b f 2 
b~ 



X = (L(l)L(2)\C\R(l)R(2)) = l 2 U 



(A8) 



+ 2d f ert(Vbf) ] 



(Al) 



-do(doTf) ctf\Vb(d Q + f) 



d 



2^(1 



±(do - /) erf[Vb(d - /)] T (do + /) erf[V&(d +/)]}, 

where erf(z) is the error function given by erf(z) = 
J Q Z e _t (ii. The detuning is given by e = tR — (l'- 



and finally there is the matrix element 



w = (R(l)L(2)\C\L/R(l)L/R(2)) 
= -j=(S(l,l)\d\S(2,0)/S(0,2)) 



l c Vbe-^I f 2 



(A9) 



d 



= { e -f)(do+/) 2 _ e ~b(d a -f) 2 



In the basis built of the orthogonalized $^ /r states we 
have then the bare tunneling 



/' = L I 2d - (do - /)erf[V6(do - /)] + 2/erf[V6/] 



VHi - 1 2 ) 

+V^{(-do + f) erf[Vb(d - /)] 
+(d + /) e r f{Vb(d Q + /)]}}, 



(A2) 



2(1 -F) ^ 
-(do + /)erf[V6(do + /)] 



(A10) 



All the matrix elements of Coulomb interaction in the , . , . n . . , .... 

, m , • \ l ■ i j 1 he matrix elements of the Coulomb interaction m this 

orthogonalized (Wanmer) basis can be expressed using , . 

four integrals involving the original ip± states (we use the asis are S lven Y 

± notation interchangeably with R/L). The confinement 

potential-related contribution W v defined in Eq. [14] in 

quartic and biquadratic models are given by the following w 

formulas respectively; 



W v 
W v 



2d 



(A3) 



(l-e- w 2)+2dg(l- erf[V&d ]) (A4) 



Apart from the direct and exchange integrals Do and 
Eo from Eqs (fI51)-(IT£ 
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Da = cVbe- bd « I Q (bd 2 ) , 
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2(1 -I 2 ) 2 
I 2 [(-3 + I 2 ) Uo~D -Eq+ Awq/1] 
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, (All) 
(A12) 
(A13) 
, (A14) 
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